Constant flux relation for aggregation models with desorption and fragmentation 
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We study mass fluxes in aggregation models where mass transfer to large scales by aggregation 
occurs alongside desorption or fragmentation. Two models are considered. (1) A system of diffusing, 
aggregating particles with influx and outflux of particles (in-out model) (2) A system of diffusing 
aggregating particles with fragmentation (chipping model). Both these models can exist in phases 
where probability distributions are power laws. In these power law phases, we argue that the two 
point correlation function should have a certain homogeneity exponent. These arguments are based 
on the exact constant flux scaling valid for simple aggregation with input. Predictions are compared 
with Monte Carlo simulations. 
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I. INTRODUCTION 

A variety of aggregation-diffusion models have been 
constructed over the years by denning simple stochastic 
rules governing the evolution of a set of particles on a 
discrete lattice. Some of the interest in models of this 
type comes from the fact that they can be considered 
as minimal models of physical aggregation diffusion sys- 
tems which provide a theoretical framework in which to 
study such phenomena as flocculation and gelation in 
aerosols and emulsions [1]. In addition, certain aggre- 
gation models have been shown to be related to models 
describing seemingly unrelated things such as the geom- 
etry of river networks 2], the distribution of forces in 
granular media Q or the directed abelian sandpile model 
4]. In analysing such models, most of the theoretical ef- 
fort has focussed on determining the average mass den- 
sity, P(m,t). In situations where a mean- field descrip- 
tion is appropriate, this can be calculated by solving a 
Smoluchowski-like kinetic equation [5[. However, many 
of the applications are in low spatial dimensions where 
mean- field theory is typically inapplicable [f| 0, Hi- In 
such situations, diffusive fluctuations dominate, render- 
ing the determination of P(m, t) quite non-trivial. Never- 
theless, much progress has been made for specific models. 

Despite the fact that fluctuations lead to non-trivial 
statistics for the mass distribution, very little is known, 
even numerically, about higher order correlation func- 
tions 0, [nj. This is the issue which we would like to 
address here. We consider a subclass of aggregation- 
diffusion models, specifically those which include depo- 
sition of monomers. Such models reach a statistically 
stationary state where the input of small masses is bal- 
anced by the depletion of small masses to generate larger 
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ones by aggregation. In recent work [ll|, [ijj we argued 
that it is useful to think of this stationary state as anal- 
ogous to the stationary state of a turbulent system with 
the mass flux playing the role of the energy flux. One 
of the most important fruits of this analogy is the reali- 
sation that aggregation models having a stationary state 
with constant mass flux must satisfy an analogue of Kol- 
mogorov's 4/5 Law [lj, [lj]. This constraint, which we 
call the Constant Flux Relation (CFR), fixes exactly the 
scaling of a special correlation function of the mass dis- 
tribution, namely the one which carries the mass flux. 
The power of this result is that it does not require any 
mean-field assumptions and holds equally well in the fluc- 
tuation dominated regime. In fact, it determines exactly 
the scaling of the flux-carrying correlation function for 
a broad class of homogeneous aggregation kernels even 
if the scaling of P(m, t) itself is not known. While the 
determination of a single correlation function is a modest 
step when faced with the problem of determining the full 
statistics, it can nevertheless be a powerful marker. For 
example, knowledge of the CFR exponent allowed us to 
give a relatively simple proof of the multifractality of the 
mass distribution for constant kernel aggregation in one 
dimension [15j . 

In our earlier paper [l2j , we derived CFR for a simple 
aggregation-diffusion model which we referred to as the 
Mass Model (MM). It was a model of diffusing particles 
undergoing mass conserving aggregation along with input 
of particles of small mass. The cascade was in the mass 
space with driving at small mass scales, dissipation at 
infinity, conserved quantity being mass and aggregation 
the process transferring mass to larger mass scales. We 
also derived analogous results for a collection of other 
models, aiming to stress the ubiquity of the approach 
and the central role played by conservation laws. In this 
paper we restrict ourselves to the implications of CFR 
for aggregation models. For clarity we will restate the 
argument for the MM. In addition, we ask whether the 
argument leading to CFR for the MM can tell us anything 
when the conservation law is less obvious, or when the 
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driving and dissipation scales are not widely separated as 
in the MM. We address this question in the context of two 
models studied earlier in relation to nonequilibrium phase 
transitions [lq . These models, which we call the in-out 
model and the chipping model, are generalisations of the 
MM. In both, the simple transfer of mass in the MM is 
disrupted all mass scales by evaporation in the former 
model and by fragmentation in the latter. In addition, 
the source of small mass scales is generated from within 
for the chipping model and not controlled from outside. 
The analysis for the MM was exact. Here we proceed 
by analogy in cases where an exact approach is yet to 
be developed, backing up our heuristic arguments with 
numerical simulations. 

The rest of the paper is organised as follows. In Sec. |TT] 
the two models are defined and known results are briefly 
reviewed. Sec. Mil describes the CFR for these models, 
and predicts the scaling behaviour of the two point joint 
probability distribution function. Sec. IIVI contains the 
results of Monte Carlo simulations. Sec. [V] contains a 
summary and conclusions. 



II. MODELS AND REVIEW 

In this section, the two models- in-out model and chip- 
ping model - are defined and relevant earlier results are 
reviewed. The in-out model describes a system of dif- 
fusing, aggregating particles with influx and outflux of 
particles. The chipping model describes a closed system 
of diffusing, aggregating particles with fragmentation. 

We define the models on a one dimensional lattice 
with periodic boundary conditions; generalisations to d- 
dimensional hyper-cubic lattice is straight forward. 



A. In-out model 

Each site i of the lattice has a non-negative integer 
mass variable m t > 0. Given a certain configuration of 
masses at time t, the system evolves in an infinitesimal 
time dt as follows. A site i is chosen at random (with 
probability dt), and then the following events can occur, 
(i) Adsorption: with probability q/(p + q + 1), unit mass 
is adsorbed at site i; thus rrii — > rrii + 1. (ii) Desorption: 
if the mass rrii is greater than zero, then with probability 
pj (p+q+l), unit mass is desorbed from site i; thus rrii —> 
rrii — 1 provided rrii > 0. (iii) Diffusion and aggregation: 
with probability l/(p + q + 1) the mass rrii moves to a 
randomly chosen nearest neighbour; thus — > and 
rrii±\ — > rrii±\ + rrii. The initial condition is chosen to 
be to be one in rrii = for all i. The model has two 
parameters, p, q. We shall refer to this model as the in- 
out model. 

Let P(m, t) denote probability that a site has mass m 
at time t. The large time limit will be denoted by P(m), 
i.e., P(m) = lim t ^oo P(m, t). When the adsorption rate 



q is increased keeping the desorption rate p fixed, the sys- 
tem undergoes a nonequilibrium phase transition across 
a critical line q c (p) from a phase in which P(jn) has an 
exponential tail to one in which it has an algebraic tail 
for large mass; i.e, 

f e -m/m* w hen q < q c (p), 
m~ Tc when q = q c (p), (1) 
m~ T when q > q c (p), 

where m* is a q dependent cut-off, and r and r c are 
exponents characterising the power law decay [l6l . [l7| . 
The three phases will be called as the exponential phase 
(q < q c ), the critical phase (q — q c ) and the growing 
phase (q > q c ). In addition, it was argued that as a 
function of the small deviation q = q — q c , and large time 
t, P(m, q, t) displays the scaling form 

PMt)~±Y(ntf,£), (2) 

in terms of three unknown exponents <j), a, t c , and the 
two variable scaling function Y. The three exponents 
were determined in all dimensions fl8j . Of interest in this 
paper is the behaviour in one dimension wherein 18] 

p(W;t) ~ ' q » qcip) > (4) 

where the scaling functions f c (x), f g (x) — > x° when x — > 
and f c (x), f g (x) — > when x — > oo. 

B. Chipping model 

Each site i of the lattice has a non-negative integer 
mass variable m, > 0. Given a certain configuration 
of masses at time t, the system evolves in infinitesimal 
time dt as follows. A site i is chosen at random (with 
probability dt), and then the following events can occur, 
(i) Chipping: with probability w/(w + 1), unit mass is 
chipped out from site i and added to a neighbour; thus 
rrii —> rrii — I and m,±i — > rrii±i + l. (ii) Diffusion and ag- 
gregation: with probability 1/(1 +w), the mass rrii moves 
to a randomly chosen nearest neighbour; thus rrii — > 
and rrii±i — > rrii±\ + m,. The initial condition is chosen 
to be to be one in which density is uniform. The model 
has two parameters: p, the mean density and w, the ratio 
of the chipping rate to the hopping rate. We shall refer 
to this model as the chipping model. 

The system undergoes a transition in the p-w plane [HI 
[n|. There is a critical line p c {w) in the p-w plane that 
separates two types of asymptotic behaviours of P(m). 
For fixed w, as p is varied across the critical value p c {w), 
the large m behaviour of P(m) was found to be, 

f e- m / m " p < Pc (w), 

P(m) ^lm-- p = p c {w), (5) 

[ m T + infinite aggregate p > p c (w). 
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Thus, the tail of the mass distribution changes from 
having an exponential decay to an algebraic decay as p 
approaches p c from below. As one increases p beyond p c , 
this asymptotic algebraic part of the critical distribution 
remains unchanged but in addition an infinite aggregate 
forms. This means that all the additional mass (p — 
p c )V (where V is the volume of the system) condenses 
onto a single site and does not disturb the background 
critical distribution. This is analogous, in spirit, to the 
condensation of a macroscopic number of bosons onto the 
single k — mode in an ideal Bose gas as the temperature 
goes below a certain critical value. 

The critical density p c (w) was found to be [|o| 



p c (w) = s/w~+l- 1, 



(0) 



in all dimensions. In addition, it was argued that the 
exponent r is super-universal and was equal to 5/2 in 
all dimensions. In particular, it was argued that in one 
dimension 



V-t)~^(^), 



(7) 



where the scaling functions f c h(x) —> x° when i->0 and 
fch(x) —>■ when x — > oo. 



III. 



CONSTANT FLUX RELATION 



In this section, we summarise the CFR argument 
specifically for aggregation models. We then apply it to 
the in-out model and the chipping model in a heuristic 
way. The heuristic steps to CFR are as follows. First 
identify the conserved quantity and the space in which 
it flows. Second, use the equation of motion to write a 
Boltzmann like continuity equation for the average den- 
sity of the conserved quantity in this space. This equa- 
tion identifies a flux-carrying correlation function, IT, and 
a nonlinear coupling, T, controlling the flow among de- 
grees of freedom. Dimensional analysis may be used to 
determine the scaling of II corresponding to constant 
flux. In Ref . [l2| . we showed how to make this dimen- 
sional argument exact. In this paper we will restrict 
ourselves to aggregation problems where the conserved 
quantity is mass, or potentially a power of mass, con- 
tenting ourselves with a dimensional derivation. 

Using the methods of Refs.pl]. l2~l|. the mass density in 
the MM is controlled by an equation with the following 
structure: 



d t (mP(m,t))= J 



Tm^! ,m 2 ^m,m 1 ,m 2 dmidm2 + ■ ■ ■ (8) 



where T„ 



is the aggregation kernel and II r( 



is the flux-carrying correlation function of interest here. 



In Refs. 11|.|15| we showed that IF, 



takes the form 



n m ,mi,m 2 — mP(rn\, m,2, i)8(m — mi — 7722) where the 
quantity P(mi,m2,t) should be interpreted physically 
as the probability of finding particles of masses mi and 



mi on adjacent lattice sites. Of course, when written out 
properly, Eq.® is much more complex, with the ". . ." 
representing additional integrals of the form shown (but 
with permuted arguments), a diffusion term, a source 
term and a noise term. The diffusion term is neglected 
since it is zero on average for a statistically homogeneous 
system. Likewise the noise, having mean zero, is ne- 
glected on average. The source term is zero for large 
masses. We can then take Eq.® to define the flux of 
mass, J(m,t), in mass space: 



d m J(m,t) 



Tm^rai ,m 2 Hm.mi ,m 2 dui\dm2 • (9) 



Given that the aggregation kernel is homogeneous of de- 
gree, /3, and assuming that the flux-carrying correlation 
function is homogeneous of degree h, it follows that the 
flux itself scales as J ~ m@ +h+3 . The essential observa- 
tion is that a stationary state of Eq.© corresponds, from 
Eq.©, to a constant flux of mass, i.e. J(m) independent 
of m. This immediately gives the CFR scaling for the 
flux-carrying correlation function, H m ,m 1 ,m 2 - 



-(3-3. 



(10) 



Hence, in the stationary state, we expect that, for large 
masses, mP(mi, m2)5(m — mi — 777,2) ~ m _ ^~ 3 or, 



P(mi,m 2 ) ~ m 13 3 . 



(11) 



Suppose, instead of mass, the conserved quantity was 
a power of the mass, m 7 . Formally, this argument 
would lead us to expect that, in the stationary state, 



m 7 P(mi, m2)5{m — m\ — 777,2) 
give 



-P- 3 . This would 



P(mi, 777,2) ~ m 



-7-/3-2 



(12) 



which will be relevant to the subsequent discussion. 

In the MM, and the other cases considered in Ref. fl2| . 
the conserved quantity was obvious. In the in-out model 
and the chipping model, it is not clear what is conserved, 
and whether there is a cascade in some quantity. In both 
cases, if we look in the mass space, then there is loss 
of mass at all mass scales either through desorption or 
through chipping. In addition, the flux in the chipping 
model is created from within through the fragmentation 
process. To apply Eq. ([TT]) or Eq. ([12")) to the two models 
being considered, we do the following. Identify I by using 
the known behaviour of P(rn,t). Let P(m,t) have the 
scaling form 



P(m,i) = —F. 
K ' m T \t s 



(13) 



where the scaling function F(x) goes to a constant for 
small x and to zero for large x. We identify I to be 
/ = 77i 7 , where 7 is fixed by the condition that (I) ~ t. 
Thus, 
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(14) 
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For both the models considered /3 = 0. In one dimen- 
sion, for the in-out model 7 = 5/2 at q = q c {p) and 7 = 1 
for q 3> q c (p). For the chipping model in one dimension 
7 = 9/2 at p = p c (w). Thus, 



P(mi,m 2 ) 
P(mi,m 2 ) 
P(m 1: m 2 ) 

IV. 



{m-im-i) 9 ^ 4 , <? = <7c(p), in-out, c? = 1(15) 
(toito 2 ) 3 , q^$>q c (p), in-out, d = 1(16) 



(mima) 



-13/4 



p = p c , chipping, d = (117) 



MONTE CARLO SIMULATIONS 



In this section, we present simulation results for the 
in-out and chipping models and compare the results with 
the predictions in Eqs. (fTS")) and (|17p. We first present 
results for the in-out model. CFR predicts that the joint 
distribution function P(rrii, m,-+i) ~ (m,-m,-+i )- 9 / 4 . In 
simulations what we will measure is the quantity 



7r(m) 



CFR then predicts that 



dm 1 P(m, m ). 



7r(m) 
7r(m) 



-7/2 

TO ' , (7 



g c (p) in-out model, d 
2 , q ^> q c (p) in-out model, d = 



(18) 



= 1,(19) 
1. (20) 



We did Monte Carlo simulations on a one dimensional 
lattice of size L — 4096. The desorption rate p is fixed 
to be one. For p = 1, q c w 0.3072. In Fig. [TJ we show 
the variation of n(rn) with to for g = 0.3072. As can 
be seen, the exponent —3.5 is a good fit, thus consistent 
with CFR. Also shown in Fig. [Tj is II(to) for q — 1.0. For 
this value of q, the system is in the growing phase, and 
CFR predicts that II(to) ~ to -2 . Again, this is borne 
out by simulations (see top curve of Fig. [I}. 

For the chipping model, CFR predicts that at p = 
Pc{w), 



7r(m) 



-11/2 



chipping model, d = 1. 



(21) 



Monte Carlo simulations were done for lattices of size 
4096 and 8192. The chipping rate w = 24.0. p was 
chosen to be p = p c = 4.0. The results are shown in 
Fig. [U The data is consistent with the CFR predictions. 



the probability distributions are power laws. The CFR 
prediction was borne out by numerical simulations. How- 
ever, an analytic approach is lacking. For this purpose, 
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FIG. 1: The variation of n(m) with m is shown for different q's 
in the in-out model. The lower straight line has an exponent 
—3.5 and corresponds to q — q c (p)- The upper straight line 
has an exponent —2.0 and corresponds to the growing phase 
q > ?c(p). 
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FIG. 2: The variation of 7r(m) with m is shown for two differ- 
ent lattice sizes in the chipping model, p = p c . The straight 
line has an exponent —5.5. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we examined the consequences of a 
constant flux on two models where there were no obvi- 
ous conserved quantities and the dissipation and driving 
scales were not widely separated. The two models con- 
sidered were examples of systems undergoing nonequilib- 
rium phase transitions. CFR holds in the phases where 



one possibly needs to work with the effective field theories 
for these models. 

For the chipping model, the fact that the two point 
correlations have a different exponents from the mean 
field answer clearly shows that the super-universality of 
P(m) is a coincidence. It is an open problem as to why 
this exponent does not change with dimension. 
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